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Abstract The thermoelectric transport properties of nanostructured devices con- 
tinue to attract attention from theorists and experimentalist alike as the spatial con- 
finement allows for a controlled approach to transport properties of correlated mat- 
ter Most of the existing work, however, focuses on thermoelectric transport in the 
linear regime despite the fact that the nonlinear conductance of correlated quantum 
dots has been studied in some detail throughout the last decade. Here, we review 
our recent work on the effect of particle-hole asymmetry on the nonlinear transport 
properties in the vicinity of the strong coupling limit of Kondo-correlated quan- 
tum dots and extend the underlying method, a renormalized superperturbation the- 
ory on the Keldysh contour, to the thermal conductance in the nonlinear regime. 
We determine the charge, energy, and heat current through the nanostructure and 
study the nonlinear transport coefficients, the entropy production, and the fate of the 
Wiedemann-Franz law in the non-thermal steady-state. Our approach is based on a 
renormalized perturbation theory in terms of dual fermions around the particle-hole 
symmetric strong-coupling limit. 
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1 Introduction 



The ability to transform energy from one form to another is of great socio-economical 
importance. Electricity plays in this context a special role as modern societies tend 
to rely on its permanent availability. Yet, the efficiency with which the energy stored 
in the chemical bonds of fossil fuels is transformed into electricity is only about 30% 
while the efficiency at which photovoltaic elements turn the energy of photons into 
electricity is, at the time of writing, at a level of about 20% in commercially available 
photovoltaic cells. The major part of the stored energy ends up as heat. Utilizing part 
of this waste heat for example via the Seebeck effect in a thermoelectric generator is 
evidently of great practical interest but, as with all heat engines, the efficiency of this 
process is ultimately limited by that of the ideal Camot cycle, T7„„„, = 1 — T^m/T^„„ 
where r„,j/7;,„, is the temperature of the cold/hot reservoir respectively. The propor- 
tionality factor between the efficiency of the thermoelectric generator and that of the 
Carnot engine depends on details of charge and heat transfer processes in the heat 
engine. A quantity of interest is in this context the dimensionless figure of merit, 

- S-aT 

ZT = , (1) 

where f = (T^u + ?;,„,) /2 is the average temperature, S is the Seebeck coefficient, a 
the electrical conductivity, and K the thermal conductivity. An increase in the figure 
of merit results in an enhanced efficiency closer to 77„„. In the limit Zf — !> oo; typical 
values for Zf are of the order of Zf w 1 . 

The electrical and thermal conductivity in linear response are defined through 

I = LnVy+Li2Vr, (2) 

Q = L2iVy+L22vr. 

where J is the charge current and Q is the heat current through the system in re- 
sponse to the applied gradients in voltage (V) and temperature {T) across the sam- 
ple. The transport coefficients L,j are evaluated at equilibrium i.e. for W = 0, 
Vr — and are not entirely independent, as Onsager's relation requires that L2 i — 
TLi 2 ll38l . Onsager's relations ensure that the entropy production remains semi- 
positive definite as required by the second law of thermodynamics and are valid 
beyond the linear response regime. 

The electrical and thermal conductivity are given in terms of L, j as 

o = Ln (3) 

L2\Ln—L22L\\ 
K= , (4) 

and the Seebeck coefficient is defined by 5 = L12/L11. The definition of k and S 
reflects that both are defined for vanishing charge current I. As the transport co- 
efficients Lij are evaluated at equilibrium, the fluctuation-dissipation theorem can 
be invoked to relate the response of the system to its equilibrium fluctuation spec- 
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trum |l5]|27l . If the applied gradients in V or T are not sufficiently small, higher order 
terms will contribute significantly to I and Q resulting in nonlinear corrections to 
the electrical and thermal conductivities that require a genuine out-of-equilibrium 
treatment. A calculation of the resulting nonlinear conductivities is possible only in 
certain limiting cases. 
The Boltzmann equation, 

5F(r,p) 



dt 



.... =^VrF(r,p) + ^VpF(r,p), (5) 
collisions at at 



e.g. is a semi-classical equation for the distribution function F(r,p) in phase space 
and requires the existence of well-defined quasi-particles. In addition, further ap- 
proximations are necessary to evaluate the collision term. A frequently employed 
approximation is the relaxation time approximation which assumes that the only 
effect of the non-equilibrium situation is to drive the system back to equilibrium. 
The characteristic rate F, in which the non-equilibrium state decays is then set by 
the relaxation time T (F ^ 1/t). In the relaxation time approximation, the collision 
term is given by 

dF{r,^) F(r,p)-fi)(r,p) 

^ = J (o) 

at collisions T 

where is the equilibrium distribution function. 

For an ordinary metal, well described by Landau's phenomeno logical Fermi liq- 
uid theory, the thermal and charge transport are intimately linked as both are due to 
the same quasi-particles. This is the content of the Wiedemann-Franz law. This law 
states that in the limit of purely elastic scattering, the ratio of K and the product of 
a and T approaches a constant, 

lim = Ln , (7) 

where Lq = n'^kjj/{3e^) is the Lorenz number (kg is Boltzmann's constant and \e\ 
is the charge quantum). It is worth stressing that in general any inelastic scattering, 
e.g. with phonons or magnons may contribute to the thermal conductivity at any 
finite T: K~ + )fph„„„„ + Km-non + . . . but at T = 0, )f = (fdcar™ in a Fermi liquid. 
As a consequence of the Wiedemann-Franz law, the figure of merit, Zf , of a metal 
at sufficiently low T is determined by the thermopower (or Seebeck coefficient) S 
which is typically small. The Seebeck coefficient 5 of a simple metal can be esti- 
mated from Mott's formula ||23l . 

One possible route to obtaining higher values of ZT in metals is in utilizing regimes 
where the Wiedemann-Franz law does not hold. In a superconductor e.g. one finds 
k/ {(JT) = but the thermopower vanishes also since the flow of charge in a super- 
conductor does not give rise to a heat current. One-dimensional metals violate the 
Wiedemann-Franz law as well 1471 . In certain intermetallic rare-earth metals that 
display quantum criticality the Wiedemann-Franz law is also violated ll45l[39l . As 
the system is quantum critical, the low-lying excitations are scale-invariant and very 
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different from those of a Fermi liquid. As a result, neither the Boltzmann equation is 
applicable to treat transport due to the absence of well-defined quasi-particles, nor 
is a linear-response treatment warranted, as no intrinsic scale is present compared 
to which the applied gradients can be considered small 1251 . It therefore is to be 
expected that these systems have a rich out-of-equilibrium behavior with interesting 
thermoelectric properties ll26l . 

A particular promising route to relatively high values of Zf has been offered 
by nanostructured devices and by superlattice structures of correlated materials B6l 
[161 |29l l40l [7] |53l. Nanostructured devices also allow for a controlled way of ad- 
dressing the nonlinear transport regime. Yet, nonlinear thermal transport properties 
have so far only received limited attention. This is largely due to the lack of reliable 
methods which allow for the accurate calculation of nonlinear transport coefficients 
in strongly correlated systems. A noteable expection is some recent work on the 
nonlinear thermal transport through a molecular junction coupled to local phonons 
based on rate equations [28l. Although it remains unclear if this approach does give 
reliable transport properties at low temperatures, the authors find strong enhance- 
ment of the nonlinear transport coefficients over their linear response counterparts. 

An enhancement of the nonlinear thermoelectric transport coefficients over their 
linear-response counterparts seems natural: relaxation processes occurring at finite 
T and at finite bias voltage do not enter the transport coefficients on equal footing 
so that the breakdown of the Wiedemann-Franz law will as functions of T and at fi- 
nite non-equilibrium drive may occur differently. As a result, the nonlinear thermal 
transport regime may indeed be key in the search for optimal efficiency of thermo- 
electric heat engines. 

Here, we focus on the electronic contribution to the thermoelectric transport 
properties of strongly correlated quantum dots. In particular, we study the behav- 
ior of the heat and charge current through a quantum dot described by the single- 
level Anderson model -to be specified below- in the nonlinear transport regime. We 
study the nonlinear transport coefficients, the entropy production and the fate of the 
Wiedemann-Franz law in the nonequilibrium steady-state. In accordance with above 
arguments, we indeed find that e.g. the nonlinear thermopower is considerably en- 
hanced above its linear-response counterpart. 

The linear response regime of the single-level Anderson model has been studied 
extensively 1221 181 1411 [Tol[35l l9l. Especially Ref. |[9l gives a complete discussion of 
the linear transport properties based on the numerical renormalization group (NRG) 
method which is known to give accurate results for quantum impurity models. The 
extension of these results to the nonlinear regime is difficult as most methods that 
are able to capture the physics of strong electron correlations, like e.g. the Bethe 
Ansatz |2J, NRG ID and Quantum Monte Carlo ||43l are at present largely confined 
to thermal equilibrium. Self-consistent diagrammatic methods like the non-crossing 
approximation and perturbative schemes can be extended to the Keldysh contour to 
treat the non-equilibrium situation. These methods are conserving, as they respect 
certain Ward identities |I3|. There is however no self-consistent method that captures 
the correct groundstate of the problem f2A\. Perturbation theory in the Coulomb re- 
pulsion U on the quantum dot is in principle possible 1491 1541 . This perturbative 
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expansion can be reorganized to deal with the strong coupling problem in terms of 
renormalized parameters llTSl . As it turns out, the extension of bare perturbation 
theory in U to the Keldysh contour suffers from an artificial non-conservation of the 
charge current away from the particle-hole (p-h) symmetric point ifTTl . We recently 
proposed a scheme on the Kedysh contour that explicitly respects charge conser- 
vation even away from p-h symmetry ll32l and that builds on the classical work of 
Yamada and others B9ll54l . on Hewson's renormalized perturbation theory to treat 
the strong coupling limit, and on Oguri's extension to the p-h symmetric Anderson 
model out of equilibrium lfT8l[36l[T9ll20l . as well as on a superperturbation theory 
scheme that utilizes dual fermions B2l [TSl . This method is discussed in detail be- 
low. 

Our main purpose here is to analyze the nonlinear thermoelectric transport proper- 
ties of quantum dots whose low-energy properties are described by a single-impurity 
Anderson model, in terms of this current conserving scheme. We demonstrate that 
it is possible to have in the nonlinear regime an enhanced Seebeck coefficient and a 
reduced Wiedemann-Franz (L/Lq) ratio as compared to their linear response coun- 
terparts. 

This chapter is organized as follows. In Section |2] we discuss the issue of cur- 
rent conservation and introduce the steady-state distribution function of the spin- 
degenerate single-level Anderson model model. Section|3]gives an introduction into 
the method of ||32l . with more details in Appendix B, and sections |4}{6] discusses 
the nonlinear electric and thermoelectric transport properties of a Kondo-correlated 
quantum dot. Appendix A introduces the nonequilibrium Green functions and the 
Dyson equation on the Keldysh contour 



2 Current Conservation and the Steady State Distribution 
Function 

We are interested in describing the transport properties of a small system, i.e. a 
system with a discrete spectrum and possibly strong (local) Coulomb repulsion, 
weakly coupled to a continuum of itinerant degrees of freedom. Despite its appar- 
ent simplicity, this class of models captures very well the low-energy properties 
of many nanostructured systems ranging from semi-conductor heterostructures to 
break-junctions and molecular devices ll30l [T1 133] 1341 . 

We will concentrate on the single-impurity Anderson model (SIAM) with one local 
spin-degenerate level at energy attached to two leads (L/R) which are modeled 
in terms of non-interacting fermions and which can be held at different chemical 
potentials (jU/, and ij.r). 
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2.1 The single-impurity Anderson model out of equilibrium 

The SIAM Hamiltonian is 

H = H,+H^ + H,i_,, (8) 

with 

Hc = E E^a^a^'AAa, (9) 

Hd = Y.Eddlda + U (^d\d\ - (^44 - i 

X=L.Rk.a 

Here, He is the Hamiltonian for electrons in a single conduction band at the metallic 
leads. Hd is the Hamiltonian for localized states in the dot, including the Coulomb 
interaction, and Hd-c is the coupling term between the dot and the leads. The Cj^xu 
are fermionic operators representing the creation (annihilation) of electrons in the 
conduction band of the left (A = L) or right (A = R) metallic lead. Localized states 
at the central region (quantum dot or molecule) are represented by the fermionic 
da operators. The coefficients Vj^x represent a scattering potential which couples 
the quasi-continuum delocalized states at the leads with the localized states at the 
central region. The density of states of the leads is given by p;^ = Lk ^ (^A .k ^ ® ) 
we will assume that Pl{co) = P^lc) = p{(>^)- For simplicity, we also assume that 
in what follows p{(o) is p-h symmetric (p(— o) = p(a))) and that p-h symmetry is 
broken only locally. For notational convenience, we introduce Ed = £d + U /2, such 
that the p-h symmetric case Ed = —U /2 is simply given by Ed = 0. 

Each lead (A ^ L/R) is assumed to be in thermal equilibrium at all times and 
hence described in terms of an equilibrium distribution function with well-defined 
temperature (Ti/Tr) and chemical potential (ixl/iir), see Fig[T] The difference in 
chemical potential (/Xz, — fiR = eV) and the temperature difference (AT = Ti — Tr) 
create a particle and energy flux through the central region. 

Several analytical results are available in the literature for the p-h symmetric case 
Ed = 0, starting with the already classical series of papers by Yamada and Yosida 
and others 1521 l49l l50l BTl l54l 1211 for the equilibrium case, and extensions to the 
non-equilibrium regime by Hershfield and WiUcins [[TtI . and by Oguri 1371 . The p-h 
asymmetric system, however, has not been studied to the same extent. 

The charge current through a nanostructured object attached to non-interacting 
leads has been derived in a series of papers. One of the earliest applications of the 
Keldysh formalism in this context is the calculation of the current through a tunnel- 
ing junction by Caroli et al. |6l . A general expression for the charge current through 
an interacting region in contact with simple (i.e. non-interacting) leads follows from 
the continuity equation describing the change in particle number in the lead ||3T1 . 
As shown by Hershfield and Wilkins ITTl . the dot obeys 
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lR = 2l ^2rR[g-+{(o)[l-fRico)]~g+-ico)fR{co)], 

4 = 2| g(_2ri)[^-+(a))[l -/z,(co)] -8+-{(0)fL{(0)], (10) 

where Ir (7^) is the charge current from the right (left) to the dot. In Eq. ( fTOl i. we 
have defined 

in=-y ^ forA=L,R (11) 

corresponding to the effective tunneling rate to the metallic leads, so that 
npi{(o)\Vi\-^ in the limit of a flat band {Vi^^ = V;^) of infinite bandwidth, where 

(to) = Y.k.a S{(0~ Ej^x ) is the density of states at the X = L,R leads. 

It has also been shown in this context that the average of both currents satisfies 
the relation 

Il + Ir ^ f dco IFlFr t w NT . . n r / m 

(12) 

whereas the difference, representing the net flux of particles at the central region, is 
given by 

Ir-Il = 2 I ^[g+-{co)E-+{co)-g-+{co)E+-{co)]=Q. (13) 
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Fig. 1 The quantum dot is comprised of a spin-degenerate local level at and a Coulomb matrix 
element U. The voltage drop across the quantum dot is set by the difference in chemical potential 
of the leads, V = {^li — and the temperature drop is given by Ti — Tn, where Ti/Tn is the 

temperature in the left/right lead respectively. We choose the zero of energy at the Fermi level of 
the conducting leads at zero bias voltage. 
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In steady-state, this difference therefore has to vanish, Ir — Il — 0. This condition is 
satisfied, provided 

g+-{co)E-+{co) -g-+{a)E+-{(o) = 0. (14) 

This relation certainly holds in equilibrium, where the different components of the 
self energy and the Green functions are linked by the Fermi distribution /o(o) = 

= 27tiAeciC0)Mc0), 

Seii^}) = 2mAec,i(o)[l -/o(a))], 

L-+{co) = [i:;,(a))-i:;,(co)]/o(«), 

E+-{co) = [i:;(a)) [1 -/o(«)]. (15) 

For an interacting system out of equilibrium, as discussed by Hershfield et al. lfTTl . 
an effective distribution function can be defined as follows 

g-+{co) = 27tiA{co)Fuico), 

g+-{(0) = 2niA{(o)[l-Fu{(o)]. (16) 

Since the Keldysh-Schwinger constraints between the self-energy components are 
still satisfied for the system out of equilibrium, then in particular we have that + 
E ^ = E'' — E" . Therefore, based on this relation it is possible to define a function 
^{(o) in the following way 

E-+{(o) = [E''{(o)-E"{(o)].^{(o), 

E+-{0}) = [E''{(o)-E"{o})][l-^{(o)]. (17) 

Substituting these definitions into Eq. (fl4] i. one finds that current conservation in 
steady-state is ensured, if 

A{o})[E'\(o)-E"{(o)]{^{(o)[l-Fu{(o)]-Fu{(o)[l-^{(o)])=Q. (18) 

This expression vanishes when ^{(o) ~ Fu{co), that is, in analogy with the equilib- 
rium situation, in steady-state the Green function components and the self-energy 
components are related by the same distribution function Fu{co). As shown in Ap- 
pendix A, this is indeed the case for the SIAM in the wide-band limit: 

,^ico) = Fu{co). (19) 

That the renormalized superperturbation theory does indeed respect Eq. il9[ and 
therefore is current conserving was shown in Ref. Il32l . 

It is instructive to notice that the distribution function for the T-matrix of the 
interacting SIAM obeys JST] 
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F fr.^ rifL{(o)+rRM(o) + iE-+{co) 

^^('^) = A-I,nErH ■ (2'^ 

Here, we have defined A =71 + 1/?. If one substitutes the relation Z^'^{(o) — 
- Fu{co) into Eq. one finds 

4-/m£^(co) ■ ^^^^ 

Solving this equation for Fu{co) leads to 

F fr.\ f ir.\ rLfLi(a)+rRfR{co) 

Fu[(0) ^ feff[(0) ^ 11+7^ ■ 

Interestingly, one can arrive at this conclusion from an alternative consideration: 
The steady-state condition Il—Ir = for the SIAM with identical density of states 
of left and right leads Pl{co) = Pr{(o) = p{(o) can be written as 



0=—J depie) [\V\lM/-8") \V\lFu{e){g'- f) 

+ \y\lfR{g'~g")-\y\lPumg''-g% 

or 

= 1 fl'ep(e){/-^«}[ri/i(e)+rfi/fi(e)-(ri + r«)iv(e)], 

As p{(£)) and g^{(£),T ^V) ~ g"{(i),T ,V) are both semi -positive functions, the steady- 
state conditions is simply Eq. (|22]) . 

Note that the distribution function for the local T-matrix of the SIAM assumes 
the particularly simple form of Eq. (l22t in the wide-band limit with identical density 
of states for the left and right lead and in the absence of an external magnetic field. 



3 Superperturbation theory on the Keldysh contour 

We recently proposed a renormalized non-equilibrium superperturbation theory, in 
terms of dual fermions on the Keldysh contour ||32l . Our primary motivation was to 
address the issue of current conservation away from p-h symmetry {E^i ^ 0), Eq. (O. 
The term superperturbation theory was introduced in Ref. [ITSl . where a quantum 
impurity coupled to a discrete bath made up of a small number of bath states was 
considered as a reference system. 

Here, the central idea is to define the interacting (f/ ^ 0) p-h symmetric (iij ~ 0) 
case as a reference system. The solution of the reference system is known explic- 
itly in terms of a regular expansion in U, respectively the renormalized interaction 
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Strength m (|52] [SO] [H] IS] |2T] [181 [371 . An expansion around this reference system 
is expected to work well, as the potential scattering term is marginally irrelevant. 

The retarded local Green function a, = — !0(t — T')([t/(j(T),c/^, (t')]} near the 
strong-coupling fixed point in the presence of p-h asymmetry within renormalized 
superperturbation theory becomes j 



= {(0 + iA-E'E 




where the renormalized parameters are given by the expressions; 

u^X+iU/iTtA), 
which represents the renormalized Coulomb interaction, 

Ed = x+iEd 



representing the p-h asymmetry, and 



, (23) 



(24) 



(25) 



(26) 



being the renormalized width of the quasiparticle resonance. The renormalization 



factor for the quasiparticle Green function is given by X++' with the spin suscepti 
bility given by the result obtained by Yamada and Yosida 1521 |49l |50l 



l + (3-:?rV4) {U/tzAY 



(27) 



The parameter 



C = 3 



(28) 



with j3 = Fi/Fr, is a convenient measure of the asymmetry in the coupling to the 
leads. In particular, for symmetric coupling, j3 = 1, one has ^ = 3/4. A detailed 
derivation of the renormalized superperturbation theory around the p-h symmetric 
SIAM is presented in Appendix B. 

The local spectral function A (©, T,V) = —(1/ 7t)ImG[, within our approach is 
given by the expression 
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+eu<, 1 



a) 



7tT_ 

X 



2C /TzTeVy 



eV 



(29) 



and is the basis for calculating charge and energy current through the quantum dot. 
Here, we have defined 



~E=Ed/A 



(30) 



as the degree of p-h asymmetry with respect to the width of the resonance level. The 
distribution function F[j{(o,T,V) of the local T-matrix within our scheme is given 
by 



12 (— )X(«) + 0(f/^), 



(31) 



where /o (w) is the 2nd derivative of the Fermi function with respect to ft). Further 



more, 



nAJ 

= r-+(ft),r,y), 

establishing that our approach is indeed current conserving 



(32) 



4 Electric conductance in the nonlinear regime 

The electric current in steady-state is calculated from the particle current defined in 
Eq. ( fT2b . le = el. The electrical conductance for finite bias voltage across the leads, 
IJ-L — t^-R — eV , is defined as 



d(el) 



dV 



(33) 

AT=Q 



Notice that the definition implies the absence of a temperature difference between 
the leads, (7^ = Tr = T). This expression is calculated from Eq. (fTzt and Eq. (|29] |. 
For the purpose of comparing with existing experimental data, it can be written in 
the form |3l 

G(T,V)-G(T,0) feVV {eV\^{kBTV feV\ 
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Fig. 2 The linear-response conductance, calculated from Eq. ( 1361 . in units of twice the quantum 
of conductance as a function of temperature for various values of the strength of the renormalized 
interaction u and p-h asymmetry e. In the presence of p-h asymmetry, the zero-temperature limit 
of G{T) will be smaller than twice the quantum of conductance in accordance with Friedel's sum 
rule. The Kondo temperature T/c here has been defined as kgTfc = TtA /4. 



-CtvEj 



eV 



(34) 



The value for the conductance at zero bias voltage and at zero temperature is 



Go ^ G(r = 0, y = 0) ^ ( ^"j [1 + (1 - afe'] 



(35) 



with renormalized parameters defined in Eqs. (1241 - 130) . It is remarkable that this 
expression satisfies Friedel's sum rule up to second order in £, u, which predicts that 
the conductance maximum should be Go = {2e-/h) [sin(7r(«,/))]^ - {2e^/h){l-e^). 
The temperature dependence of the electric conductance at zero bias voltage is given 
by 



G(r,0)=Go 



(36) 



ans shown in Fig.|2]for for different values of e and m. Here, the transport coefficient 
CT is given by the expression lf32l 
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0.4 0.5 0.6 OJ 
U 

Fig. 3 Dependence of cj on the renormalized quantities i7 and e according to Eq. (|37j. In the 
strong-coupling limit, i.e. i7 ~ 1, cj is in good approximation independent of e reflecting the 
marginal irrelevance of the potential scattering terni. 



CT 



a:^ l+2M2+[(8-5i7)t7-3]e2 



l + (l-fi)^e2 



(37) 



Eq. ( [36] l can be compared with a phenomenological formula, which is often em- 
ployed when fitting experimental data in order to obtain the characteristic low- 
energy {i.e. Kondo) scale Tk, 



G{T,0) = 



Go 



(i + (2iA-i)(r/r^)2)^ 



(38) 



Here, s is a phenomenological parameter which is typically taken to be i w 0.2 lfT3ll . 
Note that cj and therefore i is a function of the renormalized interaction strength u 
and p-h asymmetry e. The variation of ct with u and e is shown in Fig. |3] 

According to Eqs. (l36T l and ( l38T l, the numerical value of the coefficient cj- away 
from the p-h symmetric point will depend on the actual definition used for the Kondo 
scale Tk. The same applies to the remaining transport coefficients of Eq. ( l34b . which 
are given within our approach by [32] 



cv = 



1 -C + «'(1/2 + C) + [(m-3)(m - l)C-3t7(M-2)-3]£ 
1 + (1-M)^e2 



?2 
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CTV = TT^ ( 2(1 - C) + u\l + C) + 3m"(- - C) - 



— ^ C,u{u - 1 ) - 20 + AAu 



93 49 

. _m2 + 46t7^ - y + - + (20 - 40t7 + 355^ - 245^ 



M* + 2t75-5^)C 



£2) (l + (l-fi)"e 



~\1 7.2 



CVE, 



CTVEj 



(l-M)e 



I +15 J 1 + (1-M)2e2' 



2{u-lf {u{3{u-l)u~4) + 3)£^+9u 



7.2 , c^r.i 



(39) 



The analytical expression obtained in Eq. ( [34l i can be compared with the 'universal' 
equation which has been applied to analyze experimental measurements lfT4ll44l of 
electrical conductance under steady-state conditions, for semiconductor heterostruc- 
tures (quantum dots) and single-molecule devices beyond linear response 



Go-G(r,y) 

CtGq 



a 



eV 
^bTk 



-ycT 



eVT 



(40) 



Despite the apparently 'universal' form of Eq. (l40t . different experimental systems 
seem to differ in the numerical values of the coefficients a and y. In particular, 
experiments in GaAs quantum dots lfT4l reported average values of = 0.1 and 
7g = 0.5, whereas for single-molecule devices ll44l . considerably smaller values of 
as = 0.01 and 75 = 0.1 were obtained. According to Eqs.(l34l and Eq. (l40l i. it is 
clear that the T/f-independent coefficients a and 7 can be expressed in terms of the 
transport coefficients in Eq. (|39] l by 



a 



(41) 



It is clear from the analytical expressions, Eq. (139) , that the numerical values of 
these coefficients are expected to depend on specific sample features, particularly 
the degree of p-h asymmetry e, as well as on the renormalized Coulomb interac- 
tion u. It is particularly noteworthy that, in agreement with Fermi liquid theory, in 
the strongly interacting (Kondo) limit (7 — > 1, the transport coefficients in Eq. (|39t 
become independent of the degree of p-h asymmetry e. 

Our analytical results Eq. (l39i explain the numerical values obtained for the 
transport coefficients in quantum dot experiments lfT4l . where for instance the set 
of parameters u = 0.45, jS = 1 and e = 0.1 yield a = 0.1 and 7 = 0.47, in good 
agreement with Reffl?]. On the other hand, our theory cannot explain the particular 
combination of values for the transport coefficients in single-molecule experiments 
ll44l . suggesting that other mechanisms not captured by the SIAM may play a role 
in those systems, such as scattering with local phonons. 
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5 Energy transport through the quantum dot and the 
steady-state entropy production rate 

So far, we have discussed the charge transport through the quantum dot. The charge 
current is well defined even in the nonlinear regime due to charge conservation: ^ + 
V • / = 0, where p is the local charge density and / the associated charge current. In 
the present geometry, the continuity equation assumes a particularly simple form 

dr{N)+lR-lL = 0, (42) 

where dt represents the partial derivative with respect to time and (A^) represents 
the average local occupation at the dot site. Clearly, the condition for steady-state 
is Il=Ir^ I- The energy current can be introduced in an analogous manner since 
it is also related to a conserved quantity. The local energy balance at the spatially 
localized region, i.e. the quantum dot, becomes 

d,{E)=4-Jl (43) 

Here, {E) represents the average local internal energy, whereas 7^, y| are the energy 
currents flowing from the left lead to the quantum dot (L), or from the quantum dot 
to the right lead (R), respectively. 

From a similar analysis as for the particle current, and taking into account that 
the flow of each quasiparticle involves transport of an energy quanta Ti(0, we have 
that the net energy currents are given by 

■/I - ^ncoi-lFL)^. {/(o)) - -/(a))]/f(a))} . (44) 

Here, we have defined /^^j (fi^) = 2/i(s) (co) — 1 as the distribution function for each 
lead. It is convenient then to split these functions in two pieces as follows 

/'K fK I f fK 
R — Jeff + J^~fY^ * ^ Jl ) ^ 

fl = fff f + {fl - /#) , (45) 

where we defined ffjr^{(o) = 2ff,ff{(a) — 1 as the distribution function at the local 
site. Applying the identity 

-g-+{a)L+-{co)], (46) 



one obtains 
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+2 / (/.(o)) -/«(«)), 

+2 / gr.«|^/[/(«) (47) 

To check the condition for steady-state in the total energy flow, we substract both 
currents to obtain 

Jl-4 = 2 1 ^hco[8+-{co)E-+{(o)-g-+{(o)E+-{(o)]. (48) 

It is important to notice that the same condition that we invoked for steady-state in 
particle flow, i.e. 

g+-ico)E-+i(o) -g-+{co)E+-{co) = 0, (49) 

indeed will also imply energy conservation in steady-state, J§ — Je = 0. As the 
renormahzed superperturbation is respecting Eq. d49] . it is an appropriate tool to 
study the nonlinear thermoelectric transport properties in a controlled fashion. 

The steady-state energy current Je = = Je through the quantum dot is finally 
given by 

Je = = J d(0^!^h(0A{(0)[fL{03)-fR{(0)]. (50) 

For the generalization of transport coefficients of Eq. (|2) other than Ln, the 
knowledge of the nonlinear heat current is required. The notion of heat current in 
spatially extended systems away from equilibrium is still a matter of debate ll48l[T2l . 
In the present case, none of these difficulties are pertinent as the setup is easily cast 
into a hydrodynamic language without any approximations. 

In the hydrodynamic regime, where the out-of-equilibrium dynamics is only due 
to low frequency and long wavelength excitations, the system is characterized by a 
few so-called slow variables that are (away from criticality and in the absence of any 
Goldstone bosons) are determined entirely through conservation laws. The resulting 
local equilibrium allows for consistent determination of the entropy current Js via 

^+V./5=^, (51) 

where ^ > is the entropy production rate. Within the hydrodynamic approach, 
^ is decomposed into the currents /, associated with the conserved quantities; 
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N 

^^Y^X.Ji. (52) 

The currents Ji can also be expressed in terms of the generalized forces Xi: 

J'^-L^ijXj, (53) 
J 

which is nothing but Eq. (|2]l and the requirement ^ > is ensured by the Onsager 
relations ||3F| . 

In the case considered here, where a spatially confined region is attached to non- 
interaction leads in thermal equilibrium, we can obtain Eqs. ( BTl i and i53[ without 
resorting to the hydrodynamic limit (but confined to the steady-state) and hence not 
only obtain the entropy production rate but also determine the transport coefficients 
Lij in the nonlinear regime. 

As each lead is characterized by an equilibrium distribution with well-defined 
temperature {TJTk) and chemical potential {^Jiir), the corresponding entropy pro- 
duction rate {^lI3^r) vanishes. Therefore, the entropy currents from the left lead 
to the dot, and from the dot to the right lead, are given by the expressions 

TlJs = Je- IJ-lIl, 

TrJI = J§-^rIr. (54) 

In steady-state, dt{N) = and dt{E) =0 and consequently the energy and particle 
currents satisfy 

Je = Je ^ -^e, 

Il=Ir=I- (55) 
According to Eq. ( |54] |. the entropy fluxes in steady-state must therefore obey 

s Tl Tl ' 

jR^J^_t^j (56) 

so that Eq. ( fSTT l in the present case reads 

d,{S)+Jl~4^^. (57) 

Therefore, in steady-state, where explicit time-dependencies vanish, (9, (S) = 0, and 
the entropy production rate at the dot is found to be 

3^ = -Jea(^)+Ia(^^), (58) 
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where for notational convenience we have defined Axff =\i/l — Yr- Eq- (ISST i together 
with Eq. (|52] i allows to identify the generalized forces X, in the present case. From 
Eq. (|58l l it follows that even under conditions where the charge current vanishes 
(/ = 0), there will be entropy generation at the local dot site (for 71 ^ Tr), 

^ = ^JeA (^i^ > 0, (59) 

thus reflecting the existence of an intrinsic dissipation mechanism in order to sustain 
the steady-state regime. We notice that after Eq. (l54t . it is possible to define the heat 
currents 

Jq = TlJs =Je- l^di, 

4 = V§ = J§ - ^IrIr, (60) 

where Jq is identified as a heat current from the left (L) lead to the quantum dot, 
or from the quantum dot into the right lead (R). In steady-state, we have that the 
heat currents are: 

Jq = Je — P-d, 

J^=Je-^IrI- (61) 

Notice that in general Jq ^ Jq (total internal energy is conserved, not just heat). 
Moreover, under steady-state conditions (Tj, = Ir = I, J^ = J§ = Je), substitution of 
Eq. dMT l into Eq. ( |43] ) yields 

d,{E)=J^-J^ + i^E-^lR)I = 0- (62) 
The corresponding expressions for the heat currents in steady-state are 

Jq'-'^^ ^ Je - Pl(r)I = J d(0-^^^{h(o-i^E(R))M(o)[fL{(o)-fR{(o)]. (63) 



At finite voltage, when jXE— piR = eV > 0, the second term in Eq. ( |62] | repre- 
sents the macroscopic electric work to sustain the current through the voltage dif- 
ference imposed, while the first is the net flow of heat at the local site, which is 
connected with entropy production and dissipation, as previously discussed. More- 
over, Eq. dSSl i for the local entropy production in steady-state can alternatively be 
expressed as 

jR jL 

3^^-8---Q (64) 
Tr Te 

which is just stating that the local entropy production at the local dot site must be 
given by the difference between the rate of entropy gain at the right lead Jq/Tr, and 
the entropy loss at the left lead, Jq/Te- 
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6 Thermoelectric transport at finite bias voltage 



Having identified the entropy production rate, the generalized forces and the heat 
currents, we are now in a position to address the nonlinear generalizations of L12 
and L22 of Eq. (|2]i. 

Thermal conductance is experimentally measured under conditions such that the 
electric current vanishes. This leads to the fairly general definition 



KiT,V) = 



dJr. 



d{AT) 



(65) 



7=0 



which is valid regardless of the thermal gradients and bias voltages being infinites- 
imal or finite, and therefore is applicable beyond the linear response regime. In the 
previous section, we obtained expressions for the heat currents in steady-state con- 
ditions, Eq. ( |63] ), j'^^^ ^ Je ~ 1^l{r)J- I" particular, if we restrict ourselves to the 
condition of a vanishing charge current (/ = 0), we have that the heat currents sat- 
isfy 7^ |/=0 = |/=o, with 



dco--^hcoA{(o,T,V)[fL- /r], 



(66) 



given by Eq. ( |63] ). Therefore, the problem of calculating the thermal conductance 
can be stated as 



K{T,V) 



d{AT) 
subject to the condition 



dco-^hcoAico, T, V) [fL - fR 



1=0 



J 



d(0-^A{(o,T,VM-fR] 



0. 



(67) 



(68) 



It is clear that the condition of vanishing particle current Eq. ( |68] ) is fulfilled when 
the thermal gradient and the bias voltage are related. This relation is explicitly given 
by the definition of the Seebeck coefficient S{T,V), 



SiT,V) 



dV 



d{AT) 



(69) 



7=0 



which will be discussed in detail in the next section. It is more convenient to ex- 
press Eq. (l67T l via the implicit function differentiation rule. Since Eq. ( |68] | defines 
an implicit functional relation V = V{AT)\j=o, one has 



dxi/iAT,V{AT)) 




dxif 




dv 


dxjf 


dAT 


7=0 


dAT 


V 


dAT 


7=0^ 



(70) 



AT 



Applied to the thermal conductance, Eq. ( |67| | becomes 
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(^1 d(o^^Ji(oA{(o,T,V)[fL-fR]^^^. (71) 

In the linear response regime, the situation is relatively simple, since it is suf- 
ficient to expand — = —{AT/T)(odf[)/d(0+ {eV)dfo/dco, and to substitute 
A{co,T,V = 0) in the integrand of Eqs.dMIl and (|69] |. As a result, the relation be- 
tween temperature gradient and voltage is linear, and given by ^^^y = S{T,0) 
in Eq. ( |69] l. 

In the nonlinear regime, however, the relation is not trivial at all, since the non- 
linear relation between bias voltage and temperature gradient which satisfies the 
zero electric current condition is implicitly given by Eq. (1681 ). In order to obtain 
explicit analytical expressions, we will resort to a simplifying assumption: in what 
follows we assume that the thermal gradient is sufficiently small to consider only 
linear terms in 47" in the current Eq. (l68T l, but we shall keep higher order terms in 
the finite bias voltage. This is equivalent to write the following approximation for 
Eq. dMll 

I^{V/e)G{T,V)+ATLi2{T,V)=Q. (72) 
Here, we have defined the coefficient 

Ln{T,V) = -r-' / dco^Aico,T,V)JJ^, (73) 

with fo{(o) = {e^^'^'I'^BT _|_ ' j-jjg Fermi-Dirac distribution. The nonlinear electri- 
cal conductance G{T,V) was already obtained in Eq. ( |40] i. We thus solve for the 
temperature gradient in Eq. (|72i : 

'''' 

The integral in Eq. (fTTT i is evaluated using the Sommerfeld expansion up to 0{T^), 
resulting in 

,,,iTX,.r-''li,.Tf^'-^^. (75) 

3 ll + rR d(0 

This expression can be rewritten in the form 



(76) 



where we introduced 



Nonlinear thermoelectric response of quantum dots 



21 



Fig. 4 Dependence of the lin- 
ear response quantity S/{kBT) 
in units of \e\A on the renor- 
malized quantities u and e. 
The thermopower S increases 
with increasing p-h asymme- 
try and vanishes when « — >■ 1 
or e — > 0. 
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and the coefficient e is defined in Eq. dSOb . 



(77) 



6.i Thermopower 

The conditions of appHcability of the Mott formula for the Seebeck coefficient of a 
metal has been discussed by Johnson and Mahan 1231 . Horvatic and Zlatic ll22l used 
the Mott formula to calculate the thermopower in the asymmetric SIAM. Mott's for- 
mula states that the Seebeck coefficient is related to the energy-dependent scattering 
relaxation time T(e): 



i\e\ d(0 



m=0 



Horvatic and Zlatic 11221 showed that the Seebeck coefficient is given by 



(78) 



(79) 



where E,j is defined as E,i ~ E,i 



:.ci IS uenneu as ii^/ = + i;^(0)|^ q, i.e. as the renormalized position of 
the virtual bound state, determined in order to satisfy the Friedel sum rule by the 
condition Eci/ A = cot(7r(«f/)). In this equation, the factor 7, which determines the 
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Fig. 5 Nonlinear ther- 
mopower: Shown is the 
enhancement factor F de- 
fined in Eq. ([82), evaluated 
at zero temperature for sim- 
plicity, versus bias voltage. 
The behavior near the strong 
coupling limit u = 0.95 for 
various values of the p-h 
asymmetry e for symmetric 
(/3 = 1) and asymmetric cou- 
plings (j8 = 20). As a result, 
the nonlinear Seebeck coef- 
ficient becomes smaller than 
the linear response value. 
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enhancement of the thermopower, is the inverse of the quasiparticle Green function 
renormaHzation factor, 7=1 — [dEj^i/ d(o]a)=o- 

Here, instead of assuming the appHcabihty of Mott's formula, we apply our pre- 
vious analysis for the nonlinear regime to obtain an analytical expression for the 
Seebeck coefficient. The temperature gradient as a function of the bias voltage at 
vanishing electric current is obtained by substituting Eq. i65[ and Eq. (l34l i into 
Eq. (l74l i. Differentiating with respect to the bias voltage, the Seebeck coefficient 
is obtained up to C>(^^y'^), 



SiT,V) 



dAT 



dV 



-I 



/=oy 
27t^ (l-M)e 
3\e\ 1 + (1 -m)2£2 V A 



knT 



C2 ( eV 



V 4 / 



ci\a 



1-2 



C2 f eV 



-1 



(80) 



At low temperatures and zero bias voltage, the full expression Eq. dSOl l obtained 
from our theory for the thermopower can be expressed in the simplified form 



S = S{T,V = 0) = 



3H^l + e2(l-fi 



-0{T^). 



(81) 



It is interesting to compare the thermopower obtained from our superperturba- 
tion theory with the result obtained by Zlatic and Horvatic ll22l . The renormalized 
resonance width is 4 = where the spin susceptibility, according to Yamada- 
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Yosida's results, is related to the renormalization factor of the quasi-particle Green 
function f ~ X++- Therefore, in the zero bias voltage Umit {i.e. linear response 
regime), our expression Eq. (|80t for the thermopower reduces to Eq. dSTl ) which is 
equivalent to the result of Horvatic and Zlatic, with the difference that the depen- 
dence on the renormalized interaction u is made explicit. As limj7^i(l — u)e = 0. 
Our results indicate that the thermopower decays to zero for a p-h symmetric sys- 
tem, as well as in the strongly interacting (Kondo) limit, in agreement with the 
previous theory by Horvatic and Zlatic ll22l . This is demonstrated in Fig.|4] where 
S/ {ksT /eA) as a function of p-h asymmetry is shown for various values of the 
renormalized interaction strength u. 

The nonlinear thermopower, Eq. (l80t . as a function of bias voltage shows a much 
richer behavior as compared to the Seebeck coefficient at low temperature. Compar- 
ing Eqs. (ISOl l and (ISTl i. one can introduce an enhancement factor via 



such, that F{T,V = 0) = 1. In Fig.|5] the behavior of F{T = 0,V) near the strong 
coupling limit is shown for different values of the p-h and lead-dot coupling asym- 
metry. In parallel to the linear response thermopower, the value of S{T,V) remains 
small near the strong coupling limit. The behavior of F{T = Q,V) in the regime, 
where both charge fluctuations and p-h asymmetry are present, is shown in Fig.|6l 
The enhancement factor in this regime changes sign as a function of bias voltage and 
for sufficiently large lead-dot asymmetry )3 can become large in magnitude. Finally, 
the nonlinear thermopower can become large in the region where p-h asymmetry 
is present and charge fluctuation are strong (as compared to e) as demonstrated in 



FiT,V)=S{T,V)/S, 



(82) 



FigH 



mopower: The enhancement 
factor F, defined in Eq. | |82K is 
evaluated at zero temperature 
for simplicity. In the regime, 
where both charge fluctua- 
tions and p-h asymmetry are 
present, the nonlinear ther- 
mopower changes sign as a 
function of bias voltage and 
for sufficient lead-dot asym- 
metry P can become large in 
magnitude. 



Fig. 6 Nonlinear ther- 
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Fig. 7 Nonlinear ther- 
mopower: Shown is the en- 
hancement factor F defined 
in Eq. ([82) (here evaluated at 
zero temperature for simplic- 
ity) versus bias voltage for 
symmetric lead-to-dot cou- 
pling /3 = 1 and for j3 = 20 
at (7 = 0.5, reflecting the pres- 
ence of charge fluctuations. 
The enhancement factor tums 
out to be positive and large 
implying a huge enhancement 
of the thermopower in the 
nonlinear regime. 




0.2 0.3 
VIA 



0.4 0.5 



6.2 Thermal Conductance 

We shall now obtain an analytical expression for the thermal conductance, according 
to Eq. dZB- We set K{T,V) = Ki{T ,V) + KiiT ,V), corresponding to the two inte- 
gral terms in Eq. dTT)- We calculate Ki{T,V), by substituting d[fL- /r]/ d{AT) = 
— {co/T)d fo/ dco in the integrand as follows 

K,{T,V) = j ^ d(0(O^A{co) (^-^ ) = -At;h — ^A{Q,VJ), (83) 

where we have used the Sommerfeld expansion for the Fermi function to evaluate 
the integral. Let us now consider the contribution K2{T,V) arising from the second 
integral expression in Eq. (iTll . 



Fig. 8 Thermal conductance 
K{T,V — 0) versus temper- 
ature in units of twice the 
quantum of thermal con- 
ductance at p-h symmetry 
£ = for various values of 
the renormalized interaction 
strength u. K{T,0) shows 
already at the linear response 
level a strong dependence 
on the lead-to-dot coupling 
asymmetry j8 {kgTi^ = kA/4). 
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K2{T,V)^S{T,V) — 



dco--^hcoA{(o,T,V)[fL-fR] 

^L + iR /AT 



(84) 



Consistent with the order of approximation of Eq. (iTOt . we expand the difference of 
the distribution functions of the leads up to first order in A T, but up to third order in 
the voltage gradient, 



AT dfo dfo l-P{eV)^d^fo 

T dco dco l+p 2 dco^ 



Substituting this expansion into Eq 

t2 



6 d(o^' 
, we obtain the following expression 
3, ...,l-j3 52A(0) 



-I' 



(85) 



where we used the Sommerfeld expansion to evaluate the integrals, and the deriva- 
tives of the local spectral function are evaluated at 0) = 0, but at finite bias voltage. It 
is clear from the prefactor T^S{T,V) = 0{T^) that K2{T,V) is of ^(r^) and hence 
is beyond the order of approximation 0{T'^) as we assumed from the beginning. We 
finally obtain for the thermal conductance 



K{T,V)^2K^T-t;\l + - 



eV 



£2 



1 A A J 3 U . 



2 r TzTeV V 
3H 42 ) 



~9 
U 



m 2 



(86) 



Here, we have expressed the result in terms of the universal quantum of thermal 
conductance, Kj = n^k^T / {3h). Notice that in the linear response regime, evaluat- 
ing Eq. (1861 ) at zero bias voltage, we have 



K{T,0)^2f4- 



lai + '^w/A)^) 



1 \ 2 



(87) 



1 1 / KkeT 



From this later result, we see that as T ^ 0, we have 
K{T,0) 4C/3 

Kf^ 1+£2(1-M 



l + 'irinkBT/Ay 



(88) 
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where the factor of 2 accounts for the two independent spin "channels" of ther- 
mal conductance. In the p-hole symmetric case (e^ = 0) with symmetric contact 
couplings = 3/4), the thermal conductance at very low T is just Kt 2Kj, ac- 
counting for two universal quanta of thermal conductance. In Fig. [SI the thermal 
conductance K{T,V = 0) versus temperature is shown at the p-h symmetric point. 
In this case, the zero temperature limit only depends on the lead-to-dot coupling 
asymmetry. Fig.|9]addresses the dependence of K{T,V = 0) on p-h asymmetry as a 
function of temperature. 



6.3 Breakdown of the Wiedemann-Fmnz law 

At this level it is interesting to check the range of applicability of the Wiedemann- 
Franz law. The Wiedemann-Franz law is not expected to hold when inelastic scat- 
tering is present, which happens at finite temperature and is also expected to occur 
away from thermal equilibrium. 

We calculate the Lorenz number L by taking the ratio of Eq. ( |86] | for the thermal 
conductance, over Eq. ( [34l i for the zero-voltage electrical conductance 

K(T,V) K(T,0) 

L(r,y) = ^A^^^ (89) 
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= Lo- 



3 \A 



1 + 



Here, Lq — ^'^i/ (3e^) is the Lorenz number of a free electron gas. From Eq. 
we easily check that 

L(T) 

lim-^ = l, (90) 

which shows that the Wiedemann-Franz law is satisfied in the limit of zero tem- 
perature, independent on the strength of the interactions u and of p-h asymmetry 
e. 

At finite temperatures, inelastic scattering becomes possible and the Wiedemann- 
Franz ratio starts to deviate from the Lorenz number, as shown in Fig.[TO]for various 
values of u at the p-h symmetric point e 0. The dependence of p-h asymmetry as 
a function of u is shown for fixed temperature in Fig. [TT] 

In the introduction, we argued that dissipative processes in the nonlinear regime 
will lead to deviations from the Wiedemann-Franz law even at zero temperature. 
That this is indeed the case and that the behavior of the Wiedemann-Franz ratio 
L{T,V)/Lq away from thermal equilibrium is particularly rich, is demonstrated in 
Fig.[l2](a)-(d), where L/Lq versus bias voltage eV is shown for different regimes and 
at r = (in (a) and (b)) and atT = 0.25Tk (in (c) and (d)) for various values of lead- 




Fig. 10 Wiedemann-Franz law: The ratio L{T,V = 0)/L() approaches unity in the limit of only 
elastic scattering, i.e. at zero temperature. This is demonstrated here for various renormalized in- 
teraction strengths at p-h symmetry (e = 0). 
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U 

Fig. 11 Wiedemann-Franz law: The ratio L{T,V = 0)/Lo "is a function of the renomialized inter- 
action strength S for various values of the p-h asymmetry e at T = 0.25Tk. The effect of the p-h 
asymmetry is iirelevant in the close vicinity of the strong-coupling limit ii — ^ 1. 

to-dot coupling asymmetry. The most interesting feature is that the deviations from 
the Wiedemann-Franz law can lead to L{T,V)/Lo > 1 as well as L{T,V)/Lq < 1. 
Fig. [12] (a) shows that for symmetric lead-to-dot coupling, j3 = 1, p-h asymmetry 
only matters for very small u. Similar results are obtained for asymmetric lead-to- 
dot coupling, see Fig.[T2](b). At finite temperature and finite bias voltage the devia- 
tions from the Wiedemann-Franz law become more interesting: As seen in Fig. [12] 
(c) and (d), the strong coupling limit t7 « 1 implies a ratio L{T,V)/Lo > 0, that in- 
creases strongly with lead-to-dot asymmetry whereas values of m corresponding to 
intermediate coupling and the presence of charge fluctuations L{T,V)/L() goes from 
L{T,V)/Lq > 1 to L{T,V)/Lo < as a function of bias voltage. 



7 Conclusion 

Thermoelectric transport properties of interacting systems beyond the linear-response 
regime are largely unexplored. This is mainly due to the lack of reliable methods that 
can treat the interaction problem away from thermal equilibrium. Yet, there is good 
reason to believe that a better understanding of nonlinear thermoelectric transport 
will help in the search for better thermoelectrics with high figure of merits. 
To address the nonlinear transport properties of a strongly correlated system in a 
well-defined and traceable setting, we here considered the case of a singleimpurity 
Anderson model that is driven out of equilibrium by a finite voltage drop and a 
thermal gradient. Our approach extends a method we recently proposed for the cal- 
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Fig. 12 Behavior of the Wiedemann-Franz ratio L/Lo defined in Eq. ([89) in the nonlinear 
regime.for different temperatures and couplings strengths j8 as a function of applied bias volt- 
age V. (a) The Wiedemann-Franz ratio at zero temperature and equal lead-to-dot coupling of left 
and right lead, /3 = 1, for various renormalized interactions u versus p-h asymmetry e. The over- 
all behavior is rather similar and p-h asymmetry only matters at comparatively small interaction 
strengths. In the limit of vanishing bias voltage, the Wiedemann-Franz law, L/Lq = 1, is recovered, 
(b) similar to (a) but with (3 = 10. (c) similar to (a) but at finite temperature T = O.lSTg where 
kgTic = 7tA/4. The increase in L/Lo reflects the finite temperature behavior seen at eV = 0, see 
Fig. [To] (d) similar to (b) but at finite temperature T = 0.25r;i: . The Wiedemann-Franz ratio L/Lq 
displays a strong increase at finite temperature for asymmetric lead-do-couplings in the nonlinear 
regime. 



culation of the nonlinear conductance of molecular transistors and semiconductor 
quantum dots (Reference |32j|) to the calculation of thermoelectric transport coef- 
ficients. We reviewed this approach of Ref. |[32) focusing on the important issue of 
current conservation encoded in the distribution function of the dot density of states 
in presence of the leads. The entropy production rate in the non-thermal steady-state 
allowed us to generalize the linear response expressions of the thermoelectric trans- 
port coefficients. Explicit expressions for the (linear and nonlinear) thermal conduc- 
tance and thermopower are given and the breakdown of the Wiedemann-Franz law 
in the nonlinear regime is demonstrated. As demonstrated, the nonlinear regime of a 
quantum dot at intermediate coupling with charge fluctuations is characterized by an 
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enhanced Seebeck coefficient and a reduced Wiedemann-Franz ratio as compared 
to the linear-response quantities. 
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Appendix 

A Nonequilibrium Green functions 



t=— oo 




t= oo 

Fig. 13 Schwinger-Keldysh contour: The label — and + denote the time-ordered (— ) and anti- 
time-ordered (+) section of the contour. 



A Green function on the Schwinger-Keldysh contour, see Fig. [13] can be defined 

via 

G(r',0 = -'(7'f//(OfiJ(0), (91) 

where T is the time ordering operator along the Schwinger-Keldysh contour and the 
index of the field operators indicates that they are taken in the Heisenberg picture. 
This gives rise to four different functions, depending on whether t or t' is located on 
the time-ordered (— ) and anti-time-ordered {+) section of the contour. The nonequi- 
librium Green function can therefore be brought into the form 




where the index =F refers to the time-ordered (anti-time-ordered) path in the Keldysh 
contour The Dyson equation 

G = g + gi:G, (93) 

becomes a matrix equation for the selfenergy, where g is the bare Green function. 
The components of G are not independent: G^^ + G = G^ -I- G ^. The retarted 
and advanced Green functions G' and G" are defined as C — G — G ^ and 
G" = G — G+^. Similar relations links the corresponding components of E. 
In equilibrium, one finds 

g;+{(o) = M(o) [g:,,{(o)^g%{(o)] , 
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G+-{co) = [1-/0(0))] [G:„{co)-G%ico)] , 
E-+ico) = Mco) -i:;(to)] , 

E+-{co) = [1 - Mco)] - , (94) 

where /o = (exp[a)/r] + 1)^' is the Fermi function. 

The Dyson equation for the component is given by 

G+ = il+G''E'')g+-{l+E"G") + G''E+-G". (95) 

in terms of the retarded and advanced components of G and E. 
A distribution function can be defined via 

G+- ^ .^{G" -G'), (96) 
We also define a distribution function ^* for the self-energy of G^ : 

E+- = .^*{E" -E''). (97) 
The Dyson equation for G" (respectively G' ) implies 

E"-E''={G-^y-iG-^y, (98) 



or 



and therefore 



G" -G'' = G''{E" -E'')G", (99) 



G+ = ^{G" - G'') = G' jr(i:" - E'')G". (100) 

For a general initial state characterized by the distribution function defined 
by = /,„ {g" — g''), the first term of the right hand side of Eq. ( |95t vanishes 

{l+G'-E'-)g+-il+E"G") = Ml+G'-E'-)ig"~g'-){l+E"G") 

= M{l+G''E'')g"{l+E"G") - {l+G''E'-)g'{l +E"G")] 
= fi„[{l+G''E'')G"-G''il+E"G")] 
= fi„ (G" - G' - G'E^'G" + G'E'G") 
= fin{G" -G'' -G'iE" -E'')G") 

= fin{G"-G''-{G"-G''))^Q, (101) 



where equation ( |99] l was used in the last step. So, 



G+ = jr(G" - G' ) = G''^*(Z" - E'')G" (102) 
= G'E+-G\ (103) 
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where the second Une of the last equation is the Dyson equation for the com- 
ponent: G+ = G'L+-G". 
Therefore, 

Z+- = .^{Z" -Z'), (104) 
and by comparing equation ( |97] i with ( 1104b : 

^*{0},T,V)^ ^{(0j,V). (105) 

Note that the commutativity of ^ with components of G has been assumed. The 
derivation of equation ( |98] l also assumes that we start from the bare Green function 
g in the infinite past (f = — oo) and turn on the coupling to the leads and the Coulomb 
interaction U on the dot as the system evolves. Alternatively, one can start with g 
describing the m = system that is coupled to the leads. The steady-state properties 
and therefore relation (I105t do not depend on the prescription at f = — oo ifTTl . 



B Renormalized superperturbation theory on the Keldysh contour 

This appendix summarizes the renormalized superperturbation theory on the Keldysh 
contour of Ref. |[32l. The term superperturbation theory was used in Ref. |[T5l| and 
referred to a perturbation theory in terms of dual fermions around a fully interacting 
system solvable via e.g. exact diagonalization. In the renormalized superperturba- 
tion theory used here fy2\, the reference system is based on the work of Yamada 
and Yoshida [52 , 49, 50J for the symmetric Anderson model is used in the context 
of renormalized perturbation theory llTSI . 

We start from a coherent state representation of the action, i.e. 

CkXo\k,'^,(y) =ci,xo\k,'^,(y), 

da\o) =da\<y), 

(CT|4 =d*„{al (106) 

where Cj^xa ™d da are Grassmann numbers. The index X = L,R labels the two 
different leads. 

The non-equilibrium "partition function" for the system, in the Keldysh con- 
tour (see Fig. 1), is expressed in terms of a functional integral over time-dependent 



Grassmann fields, "^kXait) = ) \ ^nd ^>(f) = Jy, [ ■ Here, the indexes 

± refer to the time-ordered (-) and anti-time-ordered {+) path along the closed 
Keldysh contour. 

Z= f &[\if\\if]^[4>\4>]e'^^^'-^^^^^^\ (107) 
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i d, - ea 

idr - 



id, — Ed 
id, - Ed 



k,X.a 



Hit) 



Vki 
Va 



Here, 03 = 



1 

-1 



y.Jo,4>,{t)^'^ + is-i}'[i>',i>]. (108) 



is the third Pauli matrix. The Coulomb interaction terms are 



contained in the action S'^}' defined by 



iS'^'[i>\i>] = i dtul- i>Ua3 + ao)i>^ i>lia, + ao)i>^ 



1 r 
4 



(109) 



where Oo is the identity matrix. 

Since the action in Eq. ( 11091 ) is Gaussian in the M^li^it), WkXa{t) Grassmann 
fields, we integrate those in the partition function Eq. (1107b to obtain, in the 
frequency-space representation. 



In Eq. ( II lOl l, we defined the effective action as 



iS\<l''.<l'a 



(110) 



iS[H^,^ccA = iSum^^^acA^i I —Y.HcoEd^s'l'aco, (HI) 



where 



(112) 



is the effective action for a p-h symmetric {E^ = 0) and interacting {U ^ 0) system, 
and 
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for A =L,R 



(113) 



is the effective tunneling rate the metallic leads, which tends to — npi{(o)\Vi ^ 
in the limit of a flat band (Vj^i = ) of infinite bandwidth, with px{(o) = Y^k.a 5 (co — 
e^;L ) the density of states at the X = L,R lead. 

We bring to bear a super-perturbation scheme \A2'i to treat the term proportional 
to E(i in the effective action Eq. ( Ill lb . by using the p-h symmetric and interacting 
system described by the effective action Su [^'am > ^ura] in Eq. (1112b as a reference 
system. For that purpose, let us define co as the matrix Green function for the p-h 
symmetric reference system. 



go 



5(j,(U 6a, CO 

++ 

ta.m oa.m 



(114) 



Let us introduce the dual fermion (Grassmann) fields (paco 



fab 



7^ j where, as be- 

fore, the index =F refers to the time-ordered (anti-time-ordered) path on the Keldysh 
contour. We insert the identity. 



Det 



i'gamEdOigam) 

into the partition function, to obtain 



(115) 



(116) 



Here, we have defined Zq = Det [gaa>Ed^3gaa>]- We expand the linear terms in the 
dual fermion fields in the action Eq. (11171 ). and integrate over the fields 4>ao) to 
obtain the effective action 

iSiff[^L,^a<o] = if^^Y.^UGi^X'^oco + iSL[f,coJoco], (117) 

where the bare dual fermion Green function is defined by 

-1 



,/(0) 



gal + {gcycoEd^3gaco) 



-1 



= -ga,(o{ga,(0-E/a3) 'gCT,(a-(118) 



, -1 



One obtains a direct relation between the dual fermion Green function G^^, and the 
Green function for localized states in the dot, by noticing that the partition function 
can be written in two equivalent ways |f42l|, Eq. (II 17b and Eqs. dl 1 111 121 ). 
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xe' 



We define the matrix Daco 
of Eq ( 11201 ) to obtain. 



Ecid^, and take the functional derivative on both sides 



1 dZ 



1 5Zo 



Zq 5Dc 



(120) 



-1 



By noticing that Gac 



result Z, 



-I 5Zo 
SDaa> 



and the simple 



: fi^ ' 03, we obtain from Eq (1121b the relation. 



(121) 



In this equation, Ga.co is the Green function matrix for the interacting ([/ ^ 0) asym- 
metric {Ell 0) Anderson model, while <73 is the third Pauli matrix. In contrast, ga.a 
is the Green function for the interacting {U ^0) and symmetric (E^i = 0) Anderson 
model. Finally, G{, is the dual fermion matrix Green function, obtained from the 
solution of the matrix Dyson equation 



"CO 



'CT,(U^O-,(B"a,(U 



^f(0)yf p/(0) 



(122) 



Here, the bare dual fermion Green function is defined by Eq. (II 18b . The dual 
fermion selfenergy E(j a, is obtained from the renormalized four-point vertex of 
the reference system. This is vital for obtaining the correct behavior in the stong- 
coupling limit t7 — > 1 and for ensuring current conservation in this scheme. 

From our calculation, the retarded component of the self-energy is given by the 
expression 



(l-X++)co + E,i 



2C-4 fTtTeVy A ( U \ 




1) 



A I 



(123) 

2' 



Here, the parameter 



(124) 
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for j3 = Pi/rR, is a measure of the asymmetry in the coupUng to the leads. The 
retarded local {i.e. on the dot) Green function is therefore 



a-Ed + iA+Ejiil 1 




Xj ^^A 



3 



,(125) 



where the renormalized parameters are given by the expressions: u = X++^ I (^^)' 
which represents the renormalized Coulomb interaction, e = Ej/ A, representing 
the p-h asymmetry relative to the width of the resonance, and A = X++^^ being the 
renormaUzed width of the quasiparticle resonance. The renormalization factor for 
the quasiparticle Green function is given by X++^ with the spin susceptibility given 
by the perturbation theory result obtained by Yamada and Yosida Il52l|49ll50l . 



X++ = l + {3-n^/4){U/nA)\ 



(126) 



The corresponding renormalized local spectral function A(a), T,V) = —{1/ 7t)ImG''^ , 
within our renormalized superperturbation theory is given by 



A{coJ,V) 




(127) 



As the potential scattering term is marginally irrelevant, our perturbative ap- 
proach in Ed is expected to work well ll32l . Furthermore, the dual fermions have 
the appealing property that for Ed ^ I, the dual fermion Green function is simply 
proportional to the Green function of the reference system: 



r/(o) 



-ga,(u(gc7,(B-£'j '^3) 'ga,(B ~ -ga,(B forfi,/ > 1. (128) 



For u = 0, Eq. ( 1125b reduces to 



(129) 



reproducing the non-interacting limit exactly 



Nonlinear thermoelectric response of quantum dots 



37 



References 



1. Agralt, N., Yeyati, A.L., van Ruitenbeek, J.M.: Quantum properties of atomic-sized conduc- 
tors. Phys. Rep. 377, 81 (2003) 

2. Andrei, N., Furaya, K., Lowenstein, J.H.: Solution of the Kondo problem. Rev. Mod. Phys. 
55, 331 (1983) 

3. Baym, G., Kadanoff, L.P: Conservation laws and correlation functions. Phys. Rev. 124, 287 
(1961) 

4. Bulla, R., Costi, T.A., Pruscke, T.: Numerical renormalization group method for quantum 
impurity systems. Rev Mod. Phys. 80, 395450 (2008) 

5. Callen, H.B., Welton, T.A.: Irreversibility and generalized noise. Phys. Rev. 83, 34 (1951) 

6. Caroli, C, Combescot, R., Nozieres, P., Saint-James, D.: Direct calculation of the tunneling 
cuiTent. 1971 4, 916 (J. Phys. C) 

7. Chowdhury, I., Prasher, R., Lofgreen, K., Chrysler, G., Narasimhan, S., Mahajan, R., Koester, 
D., Alley, R., Venkatasubramanian, R.: On-chip cooling by superlattice-based thin-film ther- 
moelectrics. Nat. Nanotechnol. 4, 235 (2009) 

8. Costi, T.A., Hewson, A.C., Zlatic, V.: Transport coefficients of the Anderson model via the 
numerical renormalization group. J. Phys. C 6, 2519 (1994) 

9. Costi, T.A., Zlatic, V.: Thermoelectric transport through strongly correlated quantum dots. 
Phys. Rev B 81, 235,127 (2010) 

10. Dong, B., Lei, X.L.: Effect of the Kondo correlation on the thermopower in a quantum dot. 
J. Phys.:Condens. Matter 14, 11,747 (2002) 

11. Doyon, B., Andrei, N.: Universal aspects of nonequilibrium currents in a quantum dot. 
Phys. Rev B 73, 245,326 (2006) 

12. Dubi, Y., Ventra, M.D.: Colloquium: Heat flow and thermoelectricity in atomic and molecular 
junctions. Rev. Mod. Phys. 83, 131 (2011) 

13. Goldhaber-Gordon, D., Gores, J., Kastner, M.A., Shtrikman, H., Mahalu, D., Meirav, U.: From 
the Kondo regime to the mixed-valence regime in a single-electron transistor. Phys. Rev. Lett. 
81, 5225 (1998) 

14. Grobis, M., Rau, I.G., Potok, R.M., Shtrikman, H., Goldhaber-Gordon, D.: Universal scaling 
in nonequilibrium transport through a single channel Kondo dot. Phys. Rev. Lett. 100, 246,601 
(2008) 

15. Hafermann, H., Jung, C, Brenner, S., Katnelson, M.I., Rubtsov, A.N., Lichtenstein, A.I.: Su- 
peiperturbation solver for quantum impurity models. EPL 85, 27,007 (2009) 

16. Harman, T.C., Taylor, P.J., Walsh, M.P, LaForge, B.E.: Quantum dot superlattice themioelec- 
tric materials and devices. Science 297, 2229 (2002) 

17. Hershfield, S., Davies, J.H., Wilkins, J.: Resonant tunneling through an Anderson impurity, i. 
current in the symmetric model. Phys. Rev. B 46, 7046 (1992) 

18. Hewson, A.C.: Renormalized perturbation expansions and Fermi liquid theory. 
Phys. Rev. Lett. 70, 4007 (1993) 

19. Hewson, A.C., Bauer, J., Oguri, A.: Non-equilibrium differential conductance through a quan- 
tum dot in a magnetic field. J. Phys.:Condens. Matter 17, 5413 (2005) 

20. Hewson, A.C., Oguri, A., Bauer, J.: Renormalized perturbation approach to electron transport 
through quantum dots. In: J. Bonca, S. Kruchinin (eds.) Physical Properties of Nanosystems, 
Springer, Dordrecht, p. 10 (2010) 

21. Horvatic, B., Sokcevic, D., Zlatic, V.: Finite-temperature spectral density for the Anderson 
model. Phys. Rev. B 36, 675 (1987) 

22. Horvatic, B., Zlatic, V.: Perturbation calculation of the thermoelectric power in the asymmetric 
single-orbital Anderson model. Phys. Lett. A 73, 196 (1979) 

23. Jonson, M., Mahan, G.: Mott's formula for the thermopower and the Wiedemann-Franz law. 
Phys. Rev B 21, 4223 (1980) 

24. Kirchner, S., Kroha, J., Wolfle, P.: Dynamical properties of the Anderson impurity model 
within a diagrammatic pseudoparticle approach. Phys. Rev. B 70, 165,102 (2004) 



38 



Stefan Kirchner, Farzaneh Zamani, and Enrique Mufioz 



25. Kirchner, S., Si, Q.: Quantum criticality out of equilibrium: Steady state in a magnetic single- 
electron transistor. Phys. Rev. Lett. 103, 206,401 (2009) 

26. Kirchner, S., Si, Q.: On the concept of effective temperature in current-carrying quantum 
critical states. Phys. Status Solidi B 247, 631 (2010) 

27. Kubo, R.: Statistical-mechanical theory of irreversible processes. I. general theory and simple 
applications to magnetic and conduction problems. J. Phys. Soc. Jpn. 12, 570 (1957) 

28. Leijnse, M., Wegewijs, M., Flensberg, K.: Nonlinear thermoelectric properties of molecular 
juctions with vibrational coupling. Phys. Rev.B 82, 045,412 (2010) 

29. Majumdar, A.: Thermoelectricity in semiconductor nano structures. Science 303, 777 (2004) 

30. Manasreh, O.: Semiconductor Heterojunctions and Nanostructures. McGraw-Hill, New York 
(2005) 

31. Meir, Y., Wingreen, N.: Landauer formula for the current through an interacting electron re- 
gion. Phys. Rev. Lett. 68, 2512-2515 (1992) 

32. Mufioz, E., Bolech, C, Kirchner, S.: Universal Out-of-Equilibrium Transport in Kondo- 
Correlated Quantum Dots: Renormalized Dual Fermions on the Keldysh Contour 
Phys. Rev. Lett. 110, 016601 (2013) 

33. Natelson, D., Yu, L.H., Ciszek, J.W., Keane, Z.K., Tour, J.M.: Single-molecule transistors: 
Electron transfer in the solid state. Chem. Physics 324, 267 (2006) 

34. Nero, J.D., de Souza, F.M., Capaz, R.B.: Molecular electronics devices: Short review. J. Com- 
put. Theor. Nanosci. 7, 1 (2010) 

35. Nguyen, T.K.T., Kiselev, M.N., Kravtsov, V.E.: Thermoelectric transport through a quantum 
dot: Effects of asymmetry in Kondo channels. Phys. Rev. B 82, 113,306 (2010) 

36. Oguri, A.: Fermi-liquid theory for the Anderson model out of equlibrium. Phys. Rev. B 64, 
153,305 (2001) 

37. Oguri, A.: Out-of-equilibrium Anderson model at high and low bias voltage. J. Phys. Soc. Jpn 
74, 110(2005) 

38. Onsager, L.: Reciprocal relations in irreversible processes L Phys. Rev. 37, 405 (1931) 

39. Pfau, H., Hartmann, S., Stockert, U., Sun, P., Lausberg, S., Brando, M., Friedemann, S., Krell- 
ner, C, Geibel, C., Wirth, S., Kirchner, S., Abrahams, E., Si, Q., Steglich, F: Themial and 
electrical transport across a magnetic quantum critical point. Nature 484, 493 (2012) 

40. Poudel, B., Hao, Q., Ma, Y, Lan, Y, Minnich, A., Yu, B., Yan, X., Wang, D., Muto, A., 
Vashaee, D., Chen, X., Liu, J., Dresselhaus, M.S., Chen, G., Ren, Z.: High-thermoelectric 
performance of nanostructured bismuth antimony telluride bulk alloys. Science 320, 634 
(2008) 

41. Reguera, D., Platero, G., Bonilla, L.L., Rubi, J.M. (eds.): K. A. Matveev, Thermopower in 
Quantum Dots (1999). (Proceedings of XVI Sitges Conference on Statistical Mechanics, 
Sitges, Barcelona, Spain, 7-11 June 1999) 

42. Rubtsov, A.N., Katsnelson, M.I., Lichtenstein, A.I.: Dual fennion approach to nonlocal coiTe- 
lations in the Hubbard model. Phys. Rev. B 77, 033,101 (2008) 

43. Rubtsov, A.N., Savkin, V.V., Lichtenstein, A.I.: Continuous time quantum Monte Carlo 
method for fermions. Phys. Rev. B 72, 035,122 (2005) 

44. Scott, G.D., Keane, Z.K., Ciszek, J.W., Tour, J.M., Natelson, D.: Universal scaling of nonequi- 
librium transport in the Kondo regime of single molecule devices. Phys. Rev. B 79, 165,413 
(2009) 

45. Tanatar, M.A., Paglione, J., Petrovic, C, Taillefer, L.: Anisotropic violation of the 
Wiedemann-Franz law at a quantum critical point. Science 316, 1320 (2007) 

46. Venkatasubramanian, R., Siivola, E., Colpitts, T., O'Quinn, B.: Thin-film thermoelectric de- 
vices with high room-temperature figures of merit. Nature 413, 597 (2001) 

47. Wakeham, N., Bangura, A.F, Xu, X., Mercure, J.F., Greenblatt, M., Hussey, N.E.: Gross vio- 
lation of the Wiedemann-Franz law in a quasi-one-dimensional conductor. Nature Comm. 2, 
396 (2011) 

48. Wu, L.A., Segal, D.: Energy flux operator, current conservation and the formal Fourier's law. 
J. Phys. A: Math. Theor 42, 025,302 (2009) 

49. Yamada, K.: Perturbation expansion for the Anderson Hamiltonian. iv. Prog. Theor. Phys. 54, 
316(1975) 



Nonlinear thermoelectric response of quantum dots 



39 



50. Yamada, K.: Thermodynamical quantities in the Anderson Hamiltonian. Prog. Theor. Phys. 
55, 1345 (1976) 

51. Yamada, K.: Perturbation expansion for the asymmetric Anderson model. Prog. Theo. Phys. 
62, 354 (1979) 

52. Yosida, K., Yamada, K.: Perturbation expansion for the Anderson Hamiltonian. 
Prog. Theor. Phys. Suppl. 46, 244 (1970) 

53. Zhang, Y., Dresselhaus, M., Shi, Y, Ren, Z., Chen, G.: High thermoelectric figure-of-merit in 
Kondo insulator nanowires at low temperatures. Nano Lett. 11, 1166 (2011) 

54. Zlatic, v., Horvatic, B.: Series expansion for the symmetric Anderson Hamiltonian. 
Phys. Rev. B 28, 6904 (1983) 



